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We study a simple model for the depinning and driven steady state phases of a solid tuned 
across a polymorphic phase transition between ground states of triangular and square symmetry. 
The competition between the underlying structural phase transition in the pure system and the 
effects of the underlying disorder, as modified by the drive, stabilizes a variety of unusual dynamical 
phases. These include pinned states which may have dominantly triangular or square correlations, 
a plastically flowing liquid-like phase, a moving phase with hexatic correlations, flowing triangular 
and square states and a dynamic coexistence regime characterized by the complex interconversion 
of locally square and triangular regions. We locate these phases in a dynamical phase diagram 
and study them by defining and measuring appropriate order-parameters and their correlations. We 
demonstrate that the apparent power-law orientational correlations we obtain in our moving hexatic 
phase arise from circularly averaging an orientational correlation function which exhibits long-range 
order in the (longitudinal) drive direction and short-range order in the transverse direction. This 
calls previous simulation-based assignments of the driven hexatic glass into question. The inter- 
mediate coexistence regime exhibits several novel properties, including substantial enhancement in 
the current noise, an unusual power-law spectrum of current fluctuations and striking metastabil- 
ity effects. We show that this noise arises from the fluctuations of the interface separating locally 
square and triangular ordered regions by demonstrating a correlation between enhanced velocity 
fluctuations and local coordinations intermediate between the square and triangular. We demon- 
strate the breakdown of effective "shaking temperature" treatments in the coexistence regime by 
showing that such shaking temperatures are non-monotonic functions of the drive in this regime. 
Finally we discuss the relevance of these simulations to the anomalous behaviour seen in the peak 
effect regime of vortex lines in the disordered mixed phase of type-II superconductors. We propose 
that this anomalous behavior is directly linked to the behavior exhibited in our simulations in the 
dynamical coexistence regime, thus suggesting a possible solution to the problem of the origin of 
peak effect anomalies. 

PACS numbers: 74.25.Uv,74.25.Wx,61.43.-j,05.60.-k,64.70.K- 



I. INTRODUCTION 



The motion of an elastic medium across a quenched 
disordered background presents a simple paradigm for 
the understanding of several experiments [ij. These in- 
clude studies of the depinning of charge-density waves [H 
0] , transport measurements in the mixed phase of type-II 
superconductors [1] as well as measurements of the flow 
of colloidal particles across rough substrates 0, 0]. The 
issue of universality at continuous depinning transitions 
has traditionally dominated much of this literature, espe- 
cially in the CDW context However, the behaviour of 
non-universal quantities in the vicinity of the depinning 
transition is often of more interest to the experimenter sj . 
The nature of order, correlations and response within 
the moving phase are also questions which underly many 
recent investigations of the physics of non-equilibrium 
steady states, motivated in large part by the consider- 
able experimental literature on dynamical states of flux 



lines in the mixed phase of driven, disordered type-II 
superconductors [1] . 

The canonical example of a remarkable non-universal 
feature of the depinning transition is the peak effect of- 
ten seen in the mixed phase in the vicinity of the upper 
critical field i?c2.9]. The peak effect, a generic property 
of weakly disordered type-II superconductors, describes 
the non-monotonic behaviour of the critical current jc 
as a function of the temperature T or applied magnetic 
field H^. This critical current measures the (depin- 
ning) force required to induce observable motion of the 
vortex line array [l^ [llj . The peak effect is an often spec- 
tacular phenomenon, with jc rising sharply in a narrow 
region whose width is comparable to that of the zero- 
field superconducting transition Q. Investigations of the 
peak effect describe a host of unusual phenomena asso- 
ciated with this narrow regime 8] . These include "finger 
print phenomena" , slow voltage oscillations, history de- 
pendent dynamic response, enhanced low-frequency noise 
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with a 1//" spectrum and many other remarkable fea- 
tures. These are often collectively referred to as "peak 
effect anomalies" [ll [ll [ll [ll 111 [13 E E S, [111, 

Approaches to understanding the peak effect have typ- 
ically followed two distinct paths. The first views the 
peak effect as arising solely from the softening of the flux- 
lattice close to Hc2, as reflected in the vanishing of the 
shear elastic constant Cge ^ {H — Hc2)^, while pinning 
strengths soften more gradually, as {H — if^9)[25l|. As 
suggested initially by Pippard[23, softer lattices should 
be able to adapt better to random pinning [2^. \2^ . In 
the second class of theories, the peak effect is a reflec- 
tion of the underlying phase diagram of a weakly pinned, 
flux-line array in the - T pknel28L129|, [3^, IM 113, 
m, H M, m, M, m, M, 111, SSj^ it has been 
argued that such phase diagrams should generically ac- 
comodate intermediate glassy phases close to the melting 
transition [13, HO, [SSj, i35|, i37j. The peak effect is then 
proposed to be associated with abrupt changes in trans- 
port across such phase boundaries, with the anomalies 
rationalized in terms of the glas sy nature of such inter- 
mediate states m, m m, 133,111. 

A third, as yet unexplored, alternative to these ap- 
proaches which addresses the origin of the anomalies 
directly, combines the scenario of an underlying static 
phase transition with the possibility that driving such a 
system induces dynamical steady states with no static 
counterpart. With this motivation in mind, the central 
questions addressed in this paper are the following: Con- 
sider an underlying static phase transition in a pure sys- 
tem as modified or broadened by weak quenched disorder. 
How are signals of this transition manifest in dynamical 
measurements? Further, can novel dynamical states with 
no analog in either the pure or the disordered undriven 
system be obtained once the system is driven? Finally, 
could some of the remarkable phenomenology of the peak 
effect anomalies possibly originate in the properties of 
such states? 

We recently proposed a suitable model system capable 
of addressing some of these issues [13, EH- Our model 
uses interacting particles in two dimensions close to zero 
temperature. These particles form a crystal in the ab- 
sence of disorder. The interaction potential contains a 
simple two-body repulsive power-law interaction as well 
as a short-range three-body interaction. The two-body 
interaction favours a triangular lattice. The three-body 
term favours a square lattice. We tune between a square 
and a triangular ground state by varying the strength 
of a single parameter, the coefficient of the three-body 
interaction V3. We place N such interacting particles 
in a quenched disordered background, modeled numeri- 
cally in terms of a Gaussian random field with specified 
strength and two-point correlations. After finding the 
ground state of the interacting particles in the disordered 
background using simulated annealing techniques, we ap- 
ply a uniform driving force to the particles. This results, 
first, in a depinning transition and then a sequence of 
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FIG. 1: [Color Online] : Schematic Dynamical Phase Diagram 
of our model system, plotted as a function of the three-body 
interaction strength V3 and of driving force F (see text for a 
more detailed discussion) . The phases are (A) pinned triangu- 
lar (B) pinned square (C) plastic flow (D) anisotropic hexatic 
(E) flowing triangular (G) coexistence regime and (H) flow- 
ing square. The precise location of the lines separating these 
phases is disorder-dependent, although the general topology 
of the phase diagram is not sensitive to disorder. The bound- 
ary separating plastic flow from the flowing square phase is a 
very strong function of the driving force F at large V3. 



partially ordered states with varying degrees of spatial 
correlations as the force is increased. 

The advantages to such a formulation are several. 
Signatures of phase transitions in dynamical measure- 
ments are typically overwhelmed by thermal fluctuations 
for purely temperature driven transitions, thereby ob- 
scuring the very effects we wish to characterize. Our 
model of a T = transition between square and tri- 
angular phases surmounts this problem while also mir- 
roring similar structural phase transitions in vortex lat- 
tices, typically between triangular and distorted rectan- 
gular phaseSjacross which broadened peak effects are 
seen [ll, El, [Hi . 

The sequence of steady states obtained in our model as 
a function of the uniform external driving force F acting 
on the particles, and for various values of the three-body 
interaction strength U3, is summarized in the dynamical 
"phase" diagram of Fig. [1] This phase diagram extends 
a similar phase diagram proposed earlier to much larger 
values of U3 [3, [Z^ . The phase diagram shows a variety of 
phases: pinned states which may have dominantly trian- 
gular or square correlations, a plastically flowing "liquid- 
like" state, a moving anisotropic hexatic phase, flowing 
triangular and square states ordered over the size of our 
simulation cell and a dynamic coexistence regime. 

This paper presents a detailed characterization of these 
states. We define and calculate appropriate correlation 
functions for n— atic (where n = 6, mostly) oricnta- 
tional order, in addition to structure factors measuring 
the distribution of particles. We capture local bond- 
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orientational order in terms of distribution functions of 
a complex number characterizing the local orientation. 
This representation is used to understand the action of 
an external force in biasing the axes of orientational or- 
der. We demonstrate that the hexatic order we measure 
arises as an artefact of averaging an anisotropic quan- 
tity, which decays either exponentially to zero or to a 
constant value along two principal directions. In the "or- 
dered" moving phases, square and triangular, shown in 
Fig. [H orientational order appears to be established over 
length scales much large than our system sizes. 

In the coexistence regime (labeled (G) in the phase dia- 
gram of Fig. [1]), we investigate the nucleation and growth 
of ordered domains of one type (square or triangular) in 
another. Such nuclei are, in general, anisotropic, forming 
along the principal crystalline directions of each phase. 
The interface region of different domains is remarkably 
dynamic. We assign the substantial increase in noise we 
see within the coexistence regime to the unusual proper- 
ties of this interface, quantifying this proposal by linking 
measures of fluctuation magnitude to coordinations in- 
termediate between square and triangular. 

We measure several quantities on the dynamical side, 
including the basic current-force relations, current statis- 
tics and current noise at many different points in the 
phase diagram. We calculate the Koshelev-Vinokur 
"shaking temperature" [sTj , defined in detail below, to 
understand how disorder-induced fluctuations in the flow 
might provide an effective, pure-system temperature in 
terms of which phase behaviour can be discussed. This 
shaking temperature, while a monotonically decreasing 
function of the applied force over much of the F — 
phase diagram, is strikingly non-monotonic in F within 
the coexistence regime. We correlate local orientational 
and density fluctuations within the coexistence state to 
understand the origins of the anomalous noise in the co- 
existence regime. We compare the phenomenology of 
the simulations, specifically relating to the coexistence 
regime, to what is seen in experiments on peak effect 
anomalies in the mixed state. We argue that this compar- 
ison, supported by general phenomenological arguments, 
suggests strongly that there may be a generic explanation 
for peak effect anomalies. 

The outline of this paper is the following. In Section 
II, we describe the model system we use, explaining our 
methodology and discussing relevant features of the sim- 
ulations and what we calculate. In Section III, we provide 
an overview of the phase diagram of the driven system, 
discussing, in particular detail, the hexatic vortex glass 
and the coexistence phase. We discuss growth and fluc- 
tuations of one phase within another, the behaviour upon 
quenching and the origins of noise in this regime. In Sec- 
tion IV, we discuss some aspects of peak effect anoma- 
lies seen in the experiments, pointing out the close rela- 
tionship between the diversity we see in our simulations 
with the experimental data. We then conjecture that 
behaviour analogous to what we see in the coexistence 
regime may be generic to all driven disordered systems 



in the vicinity of an underlying static first-order phase 
transition in the pure limit. Finally, in our concluding 
section. Section V, we summarize our results briefly and 
suggest further lines of research. 



II. THE MODEL SYSTEM AND 
METHODOLOGY 

Our model system is two-dimensional and consists of 
particles with two and three-body interactions [3 |45| . 
The three-body interaction, parametrized through a sin- 
gle parameter 1)3, tunes the system across a square- 
triangular phase transition. The total interaction energy 
for particles confined to two-dimensions and labelled by 
their position vectors is thus 

V=l/2Y,V2{n,) + l/6 V,{n,r„rk), (1) 

where = |ry| = \rj -ri\. 

We take the two-body interaction to be of the power- 
law form 

1^2(r.,)=^2(-)^^ (2) 

while the three-body term is 

Vsin^rj^rk) = vslfij sin'^{A9ijk)f jk+perimitations]. (3) 

The function = /(r^) = (r,y - tq)^ for r,j < LSctq 
and otherwise and 9ijk is the angle between Vji and Tjk- 

The two-body interaction favors a triangular ground 
state while the three-body term favors 90° and 45° bonds 
and hence a square structure. Energy and length scales 
are set using V2 = 1 and (Jq ~ I. The zero-temperature 
phase diagram for particles interacting with this poten- 
tial has been calculated in Ref. [53]- As a function of 
the parameter ^3, which measures the strength of the 
three-body term, a discontinuous transition between a 
triangular lattice, obtained for V3 < 1.5, and a square 
lattice, obtained for W3 > 1.5, is seen. A similar potential 
was used by Stillinger and Weber in an early study of 
melting of a square solid,53]. 

Particles also interact with background quenched dis- 
order in the form of a one-body Gaussian random poten- 
tial field V(j(r) with zero mean and exponentially decay- 
ing (short-range) correlations. This potential field is de- 
fined on a fine grid following a methodology due to Chud- 
novsky and Dickman_54j , and interpolation is used to find 
the value of the potential at intermediate points^. The 
disorder variance is set to v'^ = 1 and its spatial correla- 
tion length is ^ = 0.12. Larkin length estimates [1. [2^. [27j 
yield La/ a ~ 100, with a = l/p^^^ the lattice parameter, 
somewhat larger than our system size. 
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A. Methodology 

The system evolves through standard Langevin dy- 
namics; 

= Vi, 

V, = f,'"*-7v, + r + 77,(0. (4) 

Here is the velocity, f;'"* the total interaction force, 
and rji (t) the random force acting on particle i, simulating 
thermal fluctuations at temperature T. A constant force 
F — {Fx,0} drives the system. The zero-mean thermal 
noise r]i (t) is specified by 

<77,(%-(i') >=2T-fd,jSit-t'), (5) 

with T — 0.1, well below the equilibrium melting tem- 
perature of the system. The unit of time t = j<Jq/v2, 
with 7 = 1 the viscosity. 



B. Simulation details 

Our system consists of N particles, with N between 
1600 and 10, 000, in a square box at number density 
p = 1.1. Configurations obtained through a simulated 
annealing procedure are the initial inputs to our Langevin 
simulations. This annealing procedure involves equilibra- 
tion using a NVT Monte Carlo scheme at fixed density 
and within a fixed background potential of a tunable am- 
plitude. The strength of the disorder is then increased in 
steps to the working disorder strength, with the system 
equilibrated at each step for around 10^ Monte Carlo 
steps. Varying the strength of disorder in this fashion 
enables the system to converge to its true minimum en- 
ergy state more efficiently than methods which employ a 
temperature annealing schedule. 

We evolve the system using a time step of 10~'''t. The 
external force is ramped up from a starting value of 
0, with the system maintained at upto 10® steps at each 

Given the local instantaneous particle density 

p(r,t)=Y,S{r-r,{t)), (6) 

i 

we calculate a variety of structural observables at equal 
time, such as the static structure factor ^(q) defined by 

Siq) ^^exp{-iq-rij). (7) 

where = — Vj. Delaunay triangulations yield the 
probability distributions P{n) of n = 4,5,6 and 7 coor- 
dinated particles {^n^i^) — We define order 
parameters 

7^= (P(4)-P(6))/(P(4) + P(6)), (8) 



to distinguish between square and triangular phases, 

V^A = (P(6) - P(5) - P(7))/(P(6) + P(5) + P(7)), (9) 

to distinguish between liquid (disordered) and triangular 
crystals and 

tAd = (P(4) - P(5) - P(7))/(P(4) + P(5) + P(7)), (10) 

to distinguish between liquid and square crystals. In ad- 
dition, we compute the hcxatic order parameter 

V'e.j = V'6(rj) = ^exp(-i66'y), (11) 

j 

and its correlations, defined via 

geir) = {MO)Mr)) (12) 

where 9 defines the bond angle associated with the vec- 
tor connecting neighbouring particles, as measured with 
respect to an arbitrary external axis. The second mo- 
ment of the distribution of '06. i is the bond-orientational 
susceptibility. Below, we describe alternative distribu- 
tion functions which quantify the extent to which the 
axes of the crystal align along the driving force direc- 
tion. These include "Argand plots" of the distribution 
of a complex quantity which characterizes local orienta- 
tional order. The dynamical variables we study include 
the center of mass velocity Vcrm and the particle flux and 
its statistics. The centre-of-mass velocity is defined via 

^c™ = 5^1^.(0), (13) 

i 

where the brackets (.) denote an average in steady state 
and Vi{t) is the velocity of particle i at time t. We mea- 
sure the particle flux by counting the number of particles 
which cross an imaginary line crossing the a; = axis 
in a single time step and then averaging this result over 
time. This flux, essentially a current j{x,y) integrated 
over all y at flxed x, which must be independent of x in 
steady state, has fluctuations about a constant value. We 
measure and discuss the power spectrum of these fluctu- 
ations. 

Koshelev and Vinokur (KV)[5l| have suggested that 
the combination of the drive and the disorder should yield 
an effective "shaking" temperature in the moving phase. 
Such an effective temperature manifests itself in trans- 
verse and longitudinal fluctuations of the velocity. We 
calculate the KV shaking temperature T^''^^ appropri- 
ate to the drive and transverse directions, obtaining it 
from 

- (E« - <m)'>/2A^> ^ = ^.y- (14) 

i 

We measure the variation of the T'^'s as a function both 
of force and at various points in our phase diagram. 
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FIG. 2; [Color Online]: Instantaneous snapshots of steady state configurations of the system at T — 0.1 with 113 — 6.0. The 
configurations (a) - (f) are respectively for driving forces Fx = 6, 10, 20, 22, 32 and 40 respectively. The particles are colored 
according to their coordination number n — 4 (magenta), 5 (green), 6 (blue), 7 (orange) and 8 (grey), as obtained from Delaunay 
triangulations. 



The transition to the coexistence phase is marked 
by the growth of highly anisotropic, obUquely inclined 
square nuclei, in a background of approximately trian- 
gularly coordinated solid. To understand their role in 
the nucleation kinetics, we define and compute quantities 
which measure the anisotropy of such clusters. To deter- 
mine the source of excess noise in the coexistence phase, 
we measure the probability distribution of the instanta- 
neous velocity excess above the mean, as a function of the 
local coordination as well as of the local density. As we 
show below, these measurements indicate that the sub- 
stantial portion of the noise originates from regions which 
are neither wholly square or triangular but to be found 
at the interface between such locally ordered states. 



III. PHASES AND PHASE DIAGRAM OF THE 
DRIVEN SYSTEM 

A qualitative understanding of the basic structure of 
the phase diagram of Fig. [T] can be obtained from snap- 
shots of instantaneous configurations in the steady state. 
Such snapshots are shown in Fig. ^ obtained from sim- 
ulations at r = 0.1 with V3 = 6.0. The configura- 



tions labelled (a)-(f) are for driving forces Fx{= F) = 
6, 10, 20, 22, 32 and 40 respectively. Particles are colored 
according to their coordination number n — A (magenta) , 
5 (green), 6 (blue), 7 (orange) and 8 (grey), obtained from 
Delaunay triangulations and a Voronoi analysis [56]. 

Disorder-induced inhomogeneities spawn local defects 
and dislocations (5 and 7 coordinated particles) in the 
system at low F. Boosting the external force trans- 
forms the disordered pinned phase, into a plastically flow- 
ing disordered fluid-like "plastic-flow" phase. This then 
turns into a coherently moving triangular phase at still 
larger F. 

We show the evolution of the coordination number his- 
tograms of the system for the corresponding forces in 
Fig. [31 In the disordered state, 5-fold and 7-fold coor- 
dinated particles are roughly similar in number as ex- 
pected. The number of particles with 5 and 7 fold coor- 
dination decrease with the drive (Fig. [3] (a),(b)). There 
appears to be a small intermediate regime where we see 
some clustering of these defects, as in (c) and (d). Fi- 
nally, at much larger F, the number of non-hexagonal 
coordinated particles declines abruptly (Fig. [3] (e),(f)) 
and the system freezes into a triangular latttice. Similar 
observations hold for the case of scans in F at larger V3, 
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FIG. 3: [Color Online] :The coordination number probability distribution at successive stages of freezing of the driven liquid, 
at driving forces Fx = 6(a), 10(&), 20(c), 22(d), 32(e) and 40(/) and with vs — 6. These forces correspond to the forces given in 
the configuration snapshots of Fig. [2] 




FIG. 4: Structure factor 5(q) for the plastic flow state (a), anisotropic hexatic phase(b), moving triangle (c) and moving 
square solid phases (d) at — 6.0 and Fx = 10,20,60 and 140 respectively. To obtain S'(q), 50 independent configurations 
were used. The structure in (d) reflects the presence of two misoriented square crystallites. 



with the difference that the ultimate large-force state is 
the flowing square lattice. Configurations in the "coex- 
istence regime" obtained at much larger forces {F ~ 100 
for U3 = 6) are discussed separately in the sections which 
follow. 



A. Pinned Phase: 



For small F the solid is pinned. In this regime, for 
F < Fc, there is no centre of mass motion at the longest 
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times we simulate. We see transient motion in the ini- 
tial stages of the application of the force, which then dies 
down once the system optimizes its location within the 
background of pinning sites. As F is increased across 
the depinning threshold, there are long transient time- 
scales for motion to set in[4^. These time-scales appear 
to diverge as the system approaches the transition, in 
agreement with expections concerning continuous depin- 
ning transitions. However, depinning appears to be hys- 
teretic, on the length and time scales we consider i.e. 
the reverse path, from depinned to pinned states yields 
an abrupt transition between these states [3 EE]- The 
depinned state just above the transition is inhomoge- 
neous and undergoes plastic fiowjSTl Issl . consistent 
with earlier numerical work. For larger Fx the velocity 
approaches the asymptotic behaviour vcm = F^. 

The pinned phase is a phase with short-range order 
in both translational and orientational correlation func- 
tions. Correlations typically extend to about 4-8 inter- 
particle spacings at the levels of disorder we consider. At 
non-zero noise strength, strictly speaking, we would ex- 
pect an activated creep component to the motion. This 
presumably lies beyond the time scales of our simulation, 
given the relatively low temperatures (T <C Tm) we work 
at. 



B. Plastic flow phase 

Upon increasing the force, we enter a regime of sub- 
stantial plastic flowfsvl Issl . \5^ . This regime is one in 
which particle motion is extremely inhomogeneous. Sim- 
ilar plastic flow is seen in a large number of simulations 
of par ticle motion in a random pinning background [30i 
M,M,M,M,M,M,EM- what is unusual here, 
however, is the strongly non-monotonic character of the 
plastic flow boundary, as illustrated in the phase diagram 
of Fig. [T] Note that the plastic flow phase boundary in 
the F — V3 plane is largely independent of F for small 
v^. However, once the critical value of for the square- 
triangular transition is crossed, this boundary becomes 
a strong function of F, with the plastic flow region ex- 
panding considerably before it collapses again. 

The structure factor 5(g) of the plastically moving 
phase (C) obtained in a narrow region just above the de- 
pinning transition consists of liquid- like isotropic rings, as 
shown in Fig. [3]Ja). For much larger values of U3, the dis- 
ordered, pinned phase appears to depin directly into the 
highly ordered moving square lattice phase, with no trace 
of a plastic flow regime. We have searched for tetratic 
phases, with algebraically decaying tetratic correlations 
in the vicinity of this depinning transition. However, no 
such phase is apparent in our numerics. 

In the plastic flow regime, transport properties are 
noisy, reflecting the underlying highly disordered na- 
ture of the phase, in agreement with previous work on 
such plastic flow states in s yste ms without competing 
phases US [Mi [13, [H, [H, H M [fil, [H, [H. We see no 



evidence for a two-peak structure in the velocity distri- 
bution function, unlike some previous work, provided we 
average enough. Thus no particle is stationary over the 
full time-scales of our simulations. 



C. Anisotropic Hexatic Phase 

As F is increased further, we encounter a narrow 
regime in which translational correlations are short- 
ranged while angle-averaged orientational correlations 
appear to decay as power laws. Within this phase the cir- 
cular ring in S{q) concentrates into six smeared peaks, 
as shown in Fig. [D^b). The presence of 6- fold order in 
the absence of the sharp Bragg-peaks associated with 
crystalline ordering suggests that this phase may have 
hexatic orientational order. We thus tentatively ider itify 
this phase as a driven hexatic, as shown in (D)[3^, l65l |: 
the terminology "anisotropic" is justified in what follows. 
However, to strengthen this assignment, other possible 
assignments, such as to a coexistence regime known to 
plague similar analyses of hexatics in two dimensional 
systems, must be ruled out [6^. 

We compute the correlations of the hexatic order pa- 
rameter for a very large system at varying values of F, as 
shown in Fig. O where N = 10000. This figure displays 
the evolution of geir), as F is varied across the phases 
(C) (D) {E) at fixed W3 = 6. We observe a sharp 
exponential decay of hexatic correlations in (C) and find 
that at the liquid to hexatic transition for F ~ 22, the 
decay fits the universal behavior 

56(0 - (15) 

expected at the fluid hexatic transition in two- 
dimensional non-disordered fluids. This exponent is ob- 
tained in a relatively narrow range of forces. Our system 
sizes are comparable to typical sizes employed to observe 
a metastable hexatic phase in two dimensional melting 
of pure solids. 

In Fig. El^a), we compare the distribution of iV'e,!^ 
over the full system, for the plastically flowing, mov- 
ing triangular and the intermediate hexatic-glass regimes. 
The figure illustrates that in the liquid/plastic regime 
{F — 10), the distribution is primarily governed by 
non six-fold coordinated particles. In the high driving 
force regime {F — 40), six- fold coordinated particles con- 
tribute in main. This raises the question of whether the 
intervening 'hexatic' regime we observe in terms of S'(q) 
and the correlation of the hexatic order parameter (r) 
at = 24, is truly hexatic or defined by a coexistence 
of the 'solid' and 'disordered' phases. If configurations 
resemble those at solid-fluid coexistence, the distribution 
would be expected to be a sum of the disordered, rela- 
tively ordered and interface distributions (weighted with 
their relative areas). 

To settle this issue we investigate the dependence of the 
distribution of IV'e. jP on the size of the system. The dis- 
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FIG. 5: [color Online]: The hexatic correlation function ge{r) 
for F — 6, 12, 14, 16, 22, 26 and 30, each averaged over 50 
independent configurations of a A'' = 10000 particle system 
at W3 = 6. The solid line indicates the universal behavior 
geir) ~ r~^^* at the liquid to hexatic transition. 



tribution is obtained by dividing the system into a num- 
ber of blocks and computing the distribution of iV'e.iP 
within each block for every configuration. In this way, in- 
formation over many length scales can be obtained from 
the same set of configurations. This provides us with the 
distribution of |^6,iP for various fractions of the total 
system. 

We find, as shown in Fig. El^b), that apart from fi- 
nite size effects, there is not much difference between the 
distributions. Also, the distributions for the two largest 
systems, viz. 0.8L and 0.9L {L being the size of the sim- 
ulation box), coincide to within statistical errors. In the 
case of two-phase coexistence one would expect a strong 
size dependence when the size of the blocks is compara- 
ble to the size of the phase separating clusters. These 
observations rule out two-phase coexistence, at least of 
the conventional kind, as an explanation for the long ori- 
entational correlations obtained on circular averaging. 

As we show below, a more detailed characterization 
of this phase yields the following results: Long-range or- 
der in orientation is always present in the direction of 
drive, while remaining short-range in the perpendicular 
direction. When averaged over all orientations, the re- 
sulting function appears power-law-like over 2-3 decades 
i.e. over length scales addressed in virtually all simula- 
tions so far. Thus, while we call the intervening phase as 
the "anisotropic hexatic" , we diverge sharply from pre- 
vious work in our claim that the true phase is never a 
true hexatic in the real sense, since it always has long- 
range order in the drive direction. We conclude that the 
apparent hexatic correlation arises, in fact, as an arte- 
fact of circularly averaging a very anisotropic correlation 
function, with qualitatively different decays in the longi- 
tudinal and transverse directions. 



1. Analogy to the XY model 

The power-law decay of orientational correlations 
shown in Fig. [51 taken together with the scahng of |'(/'6,iP 
suggests that the driven system in the vicinity of F = 22 
might best be described as a hexatic, with power-law cor- 
relations in the local orientation but short-range trans- 
lational order. Such a state is analogous to the low- 
temperature phase of the 2-dimensional XY model, where 
the Mermin- Wagner theorem rules out long-range order 
at any finite temperature but vortex-like excitations re- 
sponsible for the transition to the disordered phase exist 
chiefly as bound vortex-anti-vortex pairs. However, the 
presence of the drive direction introduces an anisotropy 
into the system. The consequences of this anisotropy 
have not, to the best of our knowledge, been explored 
in previous studies of putative hexatic phases in driven 
disordered sohds^. 

In our studies, we have been motivated by an analogy 
to the physics of the XY model in an applied field, specif- 
ically by the intriguing possibility that the drive might 
play a role equivalent to that of the magnetic field in the 
XY model. Renormalization group studies [63] of the 2- 
dimensional XY model in an external magnetic field with 
Hamiltonian, 

TLxY = ^ JijSiSj cos{6i ~ 6j) ~ h'^Si cos{0i) (16) 

where the external field h points along the positive x- 
axis, indicate three distinct phases. These are a linearly 
confined phase, a logarithmically confined phase and a 
free vortex phase obtained as the temperature is gradu- 
ally increased. Vortex-antivortex pairs at low tempera- 
tures are linearly confined by a string of overturned spins. 
With increasing temperature these strings participate in 
a proliferation transition, but vortices still remain con- 
fined due to a residual logarithmic attraction. As the 
temperature is further increased, the vortices overcome 
this attraction and are finally deconfined. 

Since the external field breaks rotational symmetry, 
the magnetization is non-zero in all the three phases. 
Further, the free energy is argued to be smooth; thus 
the distinction between phases can only be found in the 
structure and distribution of topological defects. The 
correlation function {cos{9i) cos{9j)) is then predicted to 
saturate to a constant value (6, say) for jr^ — rj\ — > oo. 
This constant b (x since the order parameter (cos(6'i)) 
couples linearly to the field h. On the other hand, the 
correlation function {sm{9i) sin(0j)) which corresponds to 
spin fluctuations in the direction transverse to the applied 
fleld (see Eq. (6.6)), is predicted to show an exponential 
decay in all the three phases. 

2. Quantifying Local Orientational Order 

To test this possibility, we must map local orienta- 
tional order to an appropriate XY-like two-dimensional 
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FIG. 6: [Color Online:] (a)Tlie distribution of the amplitude of the local hexatic order parameter |'(/'6,iP over the system for 
N = 10000 in the plastic {Fx = 10), hexatic-glass {Fx — 24) and moving triangular {Fx = 40) regimes, with vs = 6.0 (b) 
Scaling of the distribution of |V'6,iP deep in the hexatic regime for Fx — 26, computed over fractions 0.4,0.5,0.6,0.7,0.8 and 
0.9 of the simulation box length L, with V3 = 6.0. 
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FIG. 7: Mapping to the XY model : local hexagonal configuration around the particle i (diagram on left) may be analyzed in 
terms of the phasor ■06 (r^) represented as a spin situated at (diagram on right). 



vector from which defect structures can be extracted. 
This is done through the local hexatic order parameter, 
■06,1 = V'6(ri), defined earher. It is convenient to identify 
this local quantity (a phasor) with a 'hexatic-spin' with 
components Re[ipe.i] and /m[V'6,i]j thus mapping local 
geometrical order to a soft-spin XY degree of freedom. 

We illustrate this mapping from the local orientations 
of particles in real space to XY spins in Fig. [T] Unlike 
in the XY model, such hexatic spins are not attached to 
a fixed lattice but associated with the moving particles. 
They thus encode important information concerning the 
equal-time orientational correlations in the driven sys- 
tem. Fig. [5] shows real space configurations of the parti- 
cle system with associated hexatic spins at various values 
of F, illustrating the variety of associated spin configu- 
rations present in this mapping. 

These configuration maps enable the identification of 
topological defects in the ordering in the associated XY 
model. Note the presence of locally aligned regions as 
well as vortex-like excitations of strength 1 and 1/2. 



Qualitatively, we find that in both the disordered liquid 
phase and the anisotropic hexatic phases, vortex config- 
urations in such spin configurations appear to have little 
correlation with each other, suggesting that they may be 
either unbound or relatively weakly bound at best (see 
Fig. [5] (a) and (b)); the boxed regions of these figures as 
indicated are expanded in Fig. [8^1) and (ii) for (a) and 
in (iii) for (b). However, we have not been able to es- 
tablish a quantitative distinction between the disordered 
liquid phase and the anisotropic hexatic phase using our 
simulation data. 

Some quantification is, however, possible for larger 
forces, in the flowing triangular phase, (Fig.[8](c)), where 
defects in the ordering appear to be associated with 
strips, or strings, of defect. In this case, the defect is 
locally a solid with square symmetry. (An expanded plot 
of the boxed region in Fig. [Sfc) is shown in (iv), illus- 
trating the one-dimensional character of the defect). We 
expect that such defects should be linearly bound, with 
an energy proportional to the length of the strip. The 



10 




■ I I • r ■ • 1 32 1 1 • • I I 38 1 1 I I • • I I 58 

26 28 30 32 34 4B 50 52 54 56 30 32 34 36 38 40 42 46 48 50 52 54 56 



FIG. 8: [Color Online] : Typical orientation of the 'hexatic- spins' (see text) associated with particles in the liquid (a) , anisotropic- 
hexatic (b) and triangular (c) phases, at the driving forces = 10, 24 and 40 respectively. The boxes in (a), (b) and (c) are 
expanded in (i) and (ii) (corresponding to the configurations boxed in (a)), in (iii) (corresponding to the configuration boxed 
in (b)) and (iv), corresponding to the boxed configuration in (c), to enable an easier visualization of spin configurations. The 
colors correspond to local coordination, n = 6 (blue), 5 (green), 4 (magenta) and 7 (orange). Note the presence of extended 
defects in the ordering shown in (c) , associated with "strings" with local square symmetry. In all the three pictures only a part 
of the simulation box is shown for the sake of clarity and the configurations are obtained for = 6. 



binding energy in the latter case should be proportional 
to a non-equilibrium analog of a surface tension between 
the square and the triangular crystal. This argument is 
supported by calculations of the moment of inertia ten- 
sor of the set of particles which belong to such a defect, 
averaged over configurations which contain such defects. 
The largest eigenvalue (A>) of this tensor measures the 
length of these extended string- like defects. In Fig. [HI we 
show the probability distribution P(A>) at a few differ- 
ent values of the external drive. It is clear that within 
the intermediate force regime, the distribution is expo- 
nential, implying that the energy for these excitations 
scale linearly with their size. These string like excita- 
tions align preferentially along the crystallographic axes 
of the surrounding triangular lattice. As the force is 
increased, these defects offer nucleation sites for square 
crystals which are less anisotropic. The probability dis- 
tribution then ceases to be linear in A>, leading to the 
long tail in the data shown in Fig. [5] at F = 110 . 

Fig. [To] shows the probability distribution of the 
'hexatic-spin' phasors for forces F^; = 8 (a), 22 (b), 24 (c) 
and 40 (d). While Fig. fTUT a). obtained within the plas- 
tic flow phase, appears to have a uniform distribution of 
hexatic spins with angle, Figs.fTOFb). (c) and (d) display 
substantial non-uniformity in this distribution. Note the 
following feature of Fig. [TUTc) : The mapped spins tend 
to overwhelmingly point along the drive direction in the 
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FIG. 9: [Color Online] :Probability distribution of the largest 
eigenvalue A> in a semi-logarithmic scale for external forces 
F = 30,34,36,60 and 110. It is clear that for all but the 
largest force, the probability distribution is exponential show- 
ing that the energy of string-like excitations within the trian- 
gular phase is linear with size. For the largest force which is 
inside the coexistence region (see text), the probability distri- 
bution decays slower than an exponential. 




FIG. 10: [Color Online]: Probability distribution of the 'hexatic-spin' phasors for forces Fx — 8 (a), 22 (b), 24 (c) and 40 (d). 
The plane of the Argand Diagram within the area enclosed by the points ( — 1, —1), (1, —1), (1, 1) and ( — 1, 1) is subdivided into 
200 X 200 square boxes (of width 0.01). The sum of the probability distribution over these boxes is normalized to unity. 



flowing triangular state. This should be contrasted with 
the fact that the corresponding distribution function for 
the X-Y model in zero external field is isotropic across the 
disordered to quasi-long-range-ordered transition, peak- 
ing below it at I ^6 1 7^ 0, a value independent of the phase 
angle, in a finite system. This indicates that orientational 
order in our problem is strongly biased by the drive F, if 
F is sufficiently large. The probability distribution of the 
hexatic-spins in the complex Argand plane formed by the 
R6['0e^i] and Im[i/)6,i] axes clearly describes how the ex- 
ternal symmetry-breaking field (F^) builds up anisotropy 
in the our system, thus ordering the spin orientations. At 
low drive, the anisotropy is masked by disorder-induced 
fluctuations at the scale of our simulation box. As such 
fluctuations are increasingly suppressed at higher drive 
values, the system appears to organize into a coherently 
moving lattice structure whose principal axes are biased 
by the force. The spin-spin correlation functions in our 
model. Fig. [TTT a) and Fig. [TTJ^b) , obtained via the corre- 



lation functions CLeir) (longitudinal) and Creii") (trans- 
verse), correspond to the local quantities Re['06,i] and 
Im[?/>6,i]- They are 

CLsir) = {Re[MO)]Re[Mr)]) (17) 

CtsW = {Im[MO)]Im[Mr)]) (18) 

In all the phases (C), (D) and (E) of Fig. [H the orien- 
tational correlation function CLe{r) saturates asymptot- 
ically, as expected, to a constant. In phase (D), such 
saturation is obtained only in the drive direction. The 
transition from phase (D) to phase (E), when triangular 
translational order increases continuously with F^, ap- 
pears to be smooth. The correlation function CTe{f) 
decays exponentially in all the three phases. 

Note that our mapping for orientational order and its 
correlations in the moving state is free from any under- 
lying lattice effects, in contrast to earlier studies of the 
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FIG. 11: [Color Online]: (a) The correlation functions CTe,{r) plotted for a range of force values = 12, 14, 18, 20, 22, 24, 26, 28 
and 30 and for V3 = 6.0. Note the generic exponential decay obtained, (b) The correlation function Cheir) (see text) plotted 
for a range of force values = 12, 14, 18, 20, 22, 24, 26, 28 and 30. in Ref.[6a]. 
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FIG. 12: [Color Online:] The constant part in the angular 
correlation function along the drive direction goes quadrati- 
cally with the driving field strength above Fc as predicted by 
Fertig et aZ. [gj]- The data points (circles) are for «3 — 6.0 
and the dashed line is a guide to the eye. 



XY model in a field. Therefore, any anisotropy reflected 
in the spin distribution (see Fig. llOp or their correlation 
functions reflects an intrinsic property of the particular 
driven phase. The quadratic dependence of the asymp- 
totic saturation value of CLeir) on the applied field Fx 
within the fluid (plastic flow) regime, as predicted by this 
intuitive mapping, is shown in Fig. 1121 

The results presented here address one fundamental 
issue in the literature on orientational order in driven 
disordered states, in particular the question of whether a 
quasi-long-range ordered phase, the "hexatic" can exist. 
We show here conclusively that it cannot. What exists is 
an unusual intermediate state which possesses long-range 
orientational order in the direction singled out by the 



drive whereas orientational order decays exponentially in 
the transverse direction. 



D. Phases at Large Force: Square, Triangular and 
Coexistence 

What happens at still larger force values depends on 
the value of V3 . For small V3 , the system transits directly 
from anisotopic hexatic glass to triangular moving crys- 
tal. The Bragg peaks sharpen into sharp Bragg spots 
with six-fold symmetry. At intermediate values of V3, 
the system appears to undergo an unusual transition into 
what we term a "coexistence phase" , discussed in more 
detail below. In this phase, the system has both triangu- 
lar and square domains and interconverts between them 
over a broad distribution of time-scales. Inspection of 
conflgurations suggests an analogy to equilibrium phase 
coexistence with a large, heterogeneous distribution of 
domain sizes of triangular and square regions, although 
this is an explicitly non-equilibrium system. 

At larger values of V3, the coexistence regime appears 
bounded. However, the disordered, plastic flow regime 
is observed to expand. At these values of V3 the system 
undergoes a direct transition into the square phase. For 
much larger V3 , the phase boundary between plastic and 
square phases collapses again with increasing U3, reducing 
the extent of the plastic regime. In this large V3 regime, 
the system depins discontinuously and elastically from 
a pinned to moving square crystal with no intervening 
plastic flow phase that our numerics can resolve. The 
plastic flow regime (C), as well as that of the hexatic glass 
(D) expands at larger ^3 due to the frustration of local 
triangular translational order by three-body interactions. 
On further increasing F^ , the structure obtained depends 
on the value of : for low the final crystal is triangular 
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(E) whereas for large it is square (F). We assign these 
states through a study of the structure factor S{q), as 
well as the coordination number probability distributions 
shown in Fig. [31 observing that the slighly smeared six- 
fold coordination of the hexatic glass consolidates into 
sharp Bragg-peaks across the transition into the ordered 
states, as shown in Fig. [H 



E. The Coexistence phase 

For intermediate U3 and F, the system exhibits a re- 
markable "coexistence" regime (G) best described as a 
mosaic of dynamically fluctuating square and triangular 
regions. From a direct calculation of the structure factor 
we see the simultaneous appearance of peaks correspond- 
ing to hexagonal and square order [4^. The intensity of 
the peaks from the hexagonal and square phases are com- 
parable. 

As F is increased, clusters of 4-coordinatcd parti- 
cles grow in size. The evolution of the configuration 
with increasing force is shown in Fig. I13i which illus- 
trates the square-triangle domain mosaic characteristic 
of the coexistence regime. The 6-coordinated regions de- 
crease in size and the dynamics of the interfacial region 
- with predominantly 5-coordinated particles and with 
a few isolated 7-coordinated particles - shows enhanced 
and co-ordinated fluctuations. Real space configurations 
(Fig. [131) exhibit islands of square and triangular coor- 
dination connected by interfacial regions with predom- 
inately 5 coordinated particles. The configuration, as 
viewed in the co-moving frame, is extremely dynamic, 
with islands rapidly interconverting between square and 
triangle. This interconversion has complex temporal at- 
tributes. 

As substrate randomness is averaged out due to the 
motion of the particles, the shaking temperature of the 
system can reasonably be expected to decrease. This is 
reflected in the decrease of the width in velocity compo- 
nent distributions. The "shaking temperature" predic- 
tions of Koshelev and Vinokur, would indicate a ~ 1/d 
and ~ 1 /w^ nature of the fall in the transverse and drive 
directions respectively. Our results for T^f^ are in agree- 
ment with this prediction outside the coexistence regime. 

We find, as shown in Fig. [T31 that Tj^f^ is nearly inde- 
pendent of W3. However, within the coexistence regime, 
Tj^j behaves non-monotonically. Typically, for a partic- 
ular disorder configuration and for 5.5 < W3 < 8.5, T!ff^ 
appears to increase sharply at a well defined F^, signify- 
ing the start of coexistence. Within G, T^"^^ remains high 
but drops sharply at the upper limit of G, to continue 
to follow the interrupted KV behavior. This anomalous 
enhancement of fluctuation magnitudes provides strong 
evidence for a genuine coexistence phase, since increas- 
ing the driving force would be expected to reduce current 
noise monotonically once the system depins, as observed 
in all previous simulation work on related models [60l[63|. 
The limits of the coexistence region, though sharp for any 



typical disorder realization, vary considerably between re- 
alizations. 

Clusters within the coexistence state appear through 
a nucleation and growth mechanism, as modifled by the 
anisotropy induced by the presence of the drive. Fig. 1151 
shows that for a particular disorder realization at a force 
= 98, a nucleus of square region appears which grows 
with time. In all the disorder realizations we studied, 
the square nucleus first formed in the hexagonal phase 
is anisotropic and elongated along the transverse direc- 
tion (Fig. [T5l) . a consequence of the fact that the drive 
introduces a preferred direction into the system. 

Square clusters in the coexistence regime tend to be 
less anisotropic compared to the elongated string-like de- 
fects of square within a triangular background obtained 
at smaller force values. This can be seen from Fig. [TCI 
which plots the distribution of angles made, with respect 
to the drive direction, by the principal direction of the 
square nucleus, corresponding to the larger eigenvalue of 
the moment of inertia tensor. For small F, the distribu- 
tion peaks around — ±60°, indicating that the square 
nucleus forms preferentially at a 60° angle with respect to 
the drive. For forces within the coexistence regime, how- 
ever, the distribution of this angle is smooth and has no 
sharp peaks, concomitant with the shape of the nucleus 
becoming more isotropic. 



1. Noise in the Coexistence Regime 

The interconversion between square and triangular re- 
gions leads to complex spatio-temporal behavior. To 
demonstrate this, we compute the power spectrum of 
fluctuations of the particle current. We first obtain the 
statistics of the number of particles crossing an imagi- 
nary line parallel to the transverse (y) direction per unit 
time, per unit length of the line for a particular force 
{Fx) value. The power spectrum Sfiuxif) is the Fourier 
transform of the auto-correlation function of this time se- 
ries of particle flux, averaged over different choices of the 
position of the imaginary line and over different choices 
of the time slots of observation. Considerable statistics 
were taken to ensure that all quantities were well aver- 
aged. About 5 different choices of the position of the 
transverse imaginary line were taken, in addition to 20 
different choices of time slots. Averaging is performed 
over a large span of simulation time steps (^^ 10^). 

Our results are summarized in Fig.[T71for two values of 
Fx- For Fx = 60, when the system is in the moving trian- 
gular phase, Sfiuxif) is, to a large extent, featureless and 
flat, except for large / where effects due to the streaming 
of the entire system become important. In contrast, for 
a high force Fx = 100, when the system is within the co- 
existence phase, we obtain a 1//" regime, with a ~ 1.5 
in Sfiuxif) over about a decade. In addition the particle 
current fluctuations remarkably enhanced by 3-4 orders 
of magnitude in the coexistence phase as compared to 
the triangular regime. 
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FIG. 13: [Color Online:] Evolution of configurations showing square-triangle coexistence with increasing force {a)F = 
90, (6)95, (c)lOO, (d)105, (e)llO, and (/)115, together with the computed Delaunay mesh. The particles are colored accord- 
ing to the number of neighbors n = 4 (magenta), 5 (green), 6 (blue) and 7 (orange). Particles with coordination 5 are present 
mainly in the interfacial region while those with 7 are associated with isolated dislocations. 



2. Origins of Noise in the Coexistence Regime 

To investigate the origins of the noise in the coexistence 
regime, we define the local quantity 

a, =^sin2(4%fe), (19) 
j.k 

where the summation is over a defined region {r < — 
l.2a) surrounding a particle i. Here Oijk is the bond angle 
between particle i and j and k. We choose particles j and 
k such that they are all within a specified cutoff radial 
distance from particle i. 

Fig. UMa) displays the probability distribution of = 
ili for a state in the coexistence region. There are two 
prominent peaks. The peak for small values of corre- 
sponds to the square lattice (fi = for the ideal square 
crystal), whereas the peak at higher values of Q corre- 
sponds to the triangular structure. Intermediate values 
of are obtained for particles at the interface between 
square and triangle. Fig. [Wlh) shows a snapshot of 
the system in the coexistence region with particles color 
coded according to their value of ft. 

We next examine the fiuctuations of the velocity about 
the average value {{v — Wcm)^) as a function of the local 



coordination. In Fig. UWc) this is plotted as a function 
of fi. We see that velocity fluctuations are larger in the 
interfacial region. Further, Fig. [TST d) shows a plot of 
the local volume per particle as a function of the veloc- 
ity fluctuations, showing that the relatively lower density 
of the interfacial region (due to the presence of a large 
concentration of defects) causes it to fluctuate more than 
the rest of the solid. The rapid fluctuations of the inter- 
face also results in rapid intcrconversion between square 
and triangular coordinated particles and leads also to 
enhanced fluctuations in the coordination number. Since 
the driven solid is elastically constrained to move as a 
whole in the direction of the drive without gaps or cracks, 
this sets up strong correlations among the particle tra- 
jectories. Such correlations are ignored in the theory of 
the shaking temperature of Koshelev and Vinokur. 



IV. DYNAMICAL COEXISTENCE AND PEAK 
EFFECT ANOMALIES 

In this section, we turn to a possible application of 
this model, in the context of experiments on transport 
anomalies associated with peak effect phenomena in the 
superconducting mixed phase. To recapitulate, the peak 
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FIG. 14: [Color Online]: The shaking temperatures in the 
transverse (top) and drive (bottom) directions (TJ^ and T^^) 
for different V3 as a function of the centre of mass velocity 
(vcm), showing the quantity is almost independent of V3 ex- 
cept in and around the coexistence regime, due to effects of 
fluctuation. 



effect refers to the sharp increase in the critical cur- 
rent jc in the mixed phase of a disordered type-II su- 
perconductor close to Hc2(r). This critical current is 
thus a non-monotonic function of T (or H, depending 
on which is varied experimentally), since it decreases 
steadily from its low temperature value till the onset 
temperature To of the peak effect, from whence it in- 
creases sharply over a small temperature range to its 
maximum value, obtained at temperature Tm- As T is 
raised still further, the critical current collapses again, to 
reach zero in the normal state. The temperature interval 
[To,Tm] in which the critical current rises anomalously 
is the peak effect regime. This regime is dynamically 
anomalous, displaying: (i) large current noise amplifi- 
cation at low frequency, (ii) a 1// spectrum of current 
fluctuations, which is very non-Gaussian, (iii) a "finger- 
print" effect in which apparently random spikes in the 
differential resistivity as a function of drive are retraced 
as the drive is decreased, (iv) a history-dependent dy- 



namic response (v) a memory of direction, amplitude and 
frequency of applied currents, (vi) a strong suppression 
of ac response by a dc bias as well as a variety of other 
behaviour [11, [11 [11, 111, [M E E H ill, [S, [11 111 ■ 
A very large number of experiments probing such anoma- 
lous behavior, including all those referenced above, are 
transport-based, thus serving as probes of the dynamics 
of vortices within this narrow region of parameter space. 

Interpretations of these phenomena are largely phe- 
nomenological. One particularly infiuential proposal con- 
siders an underlying order-disorder transition "contam- 
inated" by sample surfaces or "edges". Such surfaces, 
with associated surface barriers for vortex entry, pro- 
vide an entry point for vortices driven into fiow|68|. The 
surface should provide an intrinsically more disordered 
environment for vortices than the bulk, particularly in 
fairly pure samples where jc is low. Thus, vortices might 
be expected to enter through the boundaries in a highly 
disordered state, only to anneal in the nearly pure bulk, 
when a current is applied across the sample. This spatial 
separation of disordered and ordered states and the slow 
annealing of one into the other is argued to be the cen- 
tral feature underlying the anomalous behavior seen in 
the peak effect regime. Magneto-optic imaging via Hall 
bar arrays support the surface contamination scenario. 
However, such methods do not access the dynamics of 
annealing and phase transformations directly. Much re- 
cent work appears consistent with a bulk coexistence of 
disordered and ordered phase [1^, while decoration ex- 
periments see a "multi-domain" structure in the peak ef- 
fect regime [toIItH. as proposed in Refs. [HIH, [s^l and 
accessed indirectly in Refs. [s^, [3^. For related simula- 
tions, see Ref. [723. 

The edge contamination scenario implicitly assumes 
that the underlying order-disorder transition is unaf- 
fected by the drive, serving only to provide a background 
to the annealing process. However, in a generic driven 
system, the possibility that the drive has more non-trivial 
effects must be expected. In particular, the drive may al- 
ter the very nature of the driven bulk, stabilizing dynam- 
ical states that are truly non-equilibrium in character, as 
illustrated in the simulations discussed above. 

A beautiful recent experiment (Ref. [73*1) performs a 
variant of scanning probe microscopy, using a mounted 
local Hall probe. The probe is sensitive to variations 
in the local magnetic induction averaged across a meso- 
scopic scale of around a micron. The system is tuned 
across the peak effect regime and then perturbed weakly 
through a low-amplitude ac field applied from below the 
sample. The Hall probe, placed above the sample and 
linked, though lock-in techniques to the frequency of the 
ac perturbation, records a local susceptibility, indicative 
of pinning response, as a function of space and integrated 
over the thickness of the sample. The spatial resolution 
is limited by the size of the Hall probe, typically of the 
order of a micron or so in size. 

As parameters are varied across the peak effect, these 
experiments see a remarkable coexistence between a 
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FIG. 15: [Color Online:] Dynamical nucleation of a square defect against a triangular background, shown for F = 98, U3 =6.0 
and at TimeSteps = 60 x 10^ (a) , 70 x 10^ (6) , 80 x 10^ (c) , 84 x 10^ (d) , 86 x 10^ (e) , 90 x 10^ (/) , 94 x 10=* (g) , 96 x 10^ (/i) , 10^ ( i) 



strong pinning regime and a weak pinning regime. A 
complex interface is seen between these coexisting states 
with a dynamics which is exquisitely sensitive to the 
field and the disorder, Such coexistence is also a feature 
of other phenomenological approaches to this problem, 
which address transport measurements. Locally more 
disordered regions of the sample appear to nucleate more 
stable regions of strong pinning whereas small variations 
of the applied field cause large changes in the inhomoge- 
neous pinning pattern. However, the complex geometry 
of the coexisting regimes appears largely stable if the 
temperature and field are fixed, suggesting that thermal 
fiuctuations are not dominant. While vortices entering 
from the sample boundaries do appear to contribute to 



this dynamics in no small measure, there is significant 
evidence for non-trivial dynamics in the bulk, with re- 
gions of strongly pinned phase being nucleated far from 
any boundary. Thus, these experiments point to a more 
active role for the bulk than envisaged in the bound- 
ary injection scenario. In this context, the authors of 
Ref. [zl, [tI] have specifically argued that the complex 
topology of the two-phase interface should be largely re- 
sponsible for the history dependence seen in the experi- 
ments. 

How are these remarkable observations related to the 
model we study here? We suggest that the link is the 
emergence of a self-organized, disorder- stabilized, dynam- 
ically sustained drive-induced coexistence phase seen both 
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FIG. 16: [Color Online:] The distribution of the angle made by the major principal axis (the eigenvector corresponding to the 
eigenvalue A>) of the square islands within a triangular background, computed for different values of the driving force F. The 
force values are (a) F = 52, (b) F = 60, (c) F = 92 and (d) F = 110, with vs, = 6. 



in the experiments and in simulations of our model sys- 
tem. The similarities between the two are striking: First, 
the coexistence itself. Both the experiments and our sim- 
ulations here provide incontrovertible evidence for dy- 
namical states in driven disordered systems which re- 
semble phase coexistence at equilibrum phase transitions, 
with the added complication o f sp atial inhomogeneities 
due to quenched disorder [gl, [tI^ [tI, [t^]. Reasoning 
from the experiments, the complex interconversion of 
one phase into the other and the spatially inhomoge- 
neous character of the dynamics is the hallmark of vor- 
tex dynamics within the peak effect regimefTSj. This is 
precisely the situation which obtains in the simulations. 
As pointed out in the previous section, what is unusual 
about the coexistence regime is that the strong disorder- 
induced fluctuations seen and manifest in all the dynamic 
properties we measure are obtained above the depinning 
transition, surviving even at large values of the applied 
force. 

Second, the coexistence seen in the experiments is very 
disorder-sensitive [73I Fzsj. In experiments, this is mani- 



fest in terms of the complex structure of differentially 
pinned regions in the sample, presumably reflecting a 
non-trivial pinning landscape. We see similar dynamical 
behavior within the coexistence regime, with the struc- 
ture at fixed drive influenced by the underlying micro- 
scopic disorder and sensitive to even marginal changes in 
the pinning. Thus, the nature of inhomogeneities con- 
nected to dynamical phase coexistence in this model as 
well as in the experiments appears to be dictated primar- 
ily by the underlying disorder. 

Third, the slow dynamics and long relaxation times 
seen in the experiments, reflecting the complex dynamics 
of the interface separating the regimes of different pinning 
strength, is seen in our simulations as wellfsl. [79|. In our 
simulations, if the drive is switched off and the system al- 
lowed to anneal, the nuclei of one phase within the other 
live anomalously long, suggesting that the dynamics has 
stabilized long-lived metastable states, precisely as seen 
in the experiments. In Fig. [12] we show the dynamics 
of a square nucleus within a triangular background, ex- 
tracted from a typical configuration within the coexis- 
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FIG. 17: [Color Online:] The power spectrum Sfiux{f) of cur- 
rent fluctuations in the moving triangle phase (open triangles 
Fx = 60) and in the coexistence region (red circles Fx = 100), 
logarithmically binned and plotted for a range of frequencies 
below the washboard frequency. Current fluctuations in the 
coexistence region are enhanced, also showing a 1/ f^^'^ decay 
in the low frequency range. 



tence regime, when quenched to zero drive. Far from van- 
ishing over a short time scale, the nucleus appears to be 
remarkably stable out to the longest time scales accessed 
in our simulations, providing evidence that the metasta- 
bility exhibited in the coexistence regime survives even 
after the drive is turned off. These are reflected in the 
memory experiments, in which the removal of the driving 
transport current essentially appears to freeze the system 
into a metastable state from which it only recovers upon 
a reapplication of the drive. 

Fourth and finally, the noise spectra within the coexis- 
tence regime, including the anomalously large noise and 
the power-law fall-off with a 1//" spectrum is a promi- 
nent feature of the peak effect regime. Experiments see 
a power-law falloff, with 1 < a < 2, as well as a sub- 
stantial enhancement in noise power of a few orders of 
magnitude [l^, [l^. The substantial non-Gaussian fea- 
tures in the noise, as obtained in Ref. [lB|, indicate that 
a small number of fluctuators contribute; it is tempting 
to assign these to a few large domains which fluctuate 
collectively within the coexistence regime, as suggested 
by the simulations. 

The edge contamination scenario assumes that the 
drive in the bulk anneals the disordered vortices injected 
across the boundaries into the smoothly flowing state. In 
contrast, we show that in the vicinity of an underlying 
transition, the bulk flowing state is generically inhomo- 
geneous and dynamically non-trivial, thereby questioning 
this basic assumption. We believe that the edge contam- 
ination scenario should apply, more immediately, close to 
but away from the peak effect regime, where the drive in 
the bulk acts solely to smoothen flow, precisely as seen 
in a recent experiment [80| . 

We stress that the system we simulate is far from a 



literal translation of the vortex system. It lacks the 
Ko{r/X), (with A the penetration depth), interactions of 
vortex lines, which can become fairly long-range if A is 
large. In our simulations we drive transitions between the 
two different pure system phases by varying a parameter 
V3, whereas these transitions are density or temperature 
driven in the classic scenario of the peak effect in vortex 
systems. (The underlying order-disorder transition in the 
vortex system is likely related to a melting-like transition, 
at least for small H and for temperatures close to Tc and 
not a structural transition between two crystalline phases 
as in our model [STj.) The experiments are in three dimen- 
sions, whereas our simulations are two-dimensional. Our 
system shows very little variation in the depinning force 
close to the transition, and thus no apparent peak effect. 
(It is known, however, that factoring in the temperature 
dependence of the coherence length ^ and A is essential to 
obtain a peak effect in two-dimensional simulations [6^.) 

For these, as well as other reasons, our proposal for the 
origin of peak effect anomalies should be thought of as in- 
dicating a generic scenario within which much of the ob- 
served physics finds a common explanation. Clearly more 
work, including simulations of the three-dimensional case 
capable of describing vortex entanglement, is required to 
further illuminate the relationship conjectured here. 



V. SUMMARY AND CONCLUSION 

This paper has studied the complex dynamical behav- 
ior of a two dimensional driven disordered solid which 
undergoes a square to triangular structural transition as 
a single parameter is tuned. Our interest in this model 
has several origins: The problem of dynamical probes of 
an underlying static phase transition in a weakly disor- 
dered system is interesting in itself. The possibility that 
such behavior might underlie a classic and ill-understood 
problem in the literature on type-II superconductivity, 
the problem of the origin of peak-effect anomalies, pro- 
vides further motivation for this study. 

Specifically, using the initial study of Ref. [43| and 
this paper, we have demonstrated the following: First, 
the existence of a complex phase diagram containing a 
large number of states induced purely by the drive, such 
as the anisotropic hexatic state and the coexistence state. 
Second, the demonstration of dynamically unusual be- 
haviour in the coexistence region, as measured through 
several dynamical quantities, including the observation 
of 1//" correlations in the current noise, the existence 
of highly metastable states and the observation of a 
very substantial noise enhancement. Third, the break- 
down of single-particle, effective temperature descrip- 
tions in the coexistence regime, where measures of a flow 
and disorder induced temperature indicate substantially 
non-monotonic variation within the coexistence regime. 
Fourth, the demonstration that the putative driven hex- 
atic glass phase cannot have algebraically decaying ori- 
entational correlations, since correlations along the drive 



19 




(b) 



(d) 




FIG. 18; [Color Online:] (a) Probability distribution of the bond angle order parameter f2 (see text) in the coexistence regime. 
The peak at low value of f2 corresponds to the square phase. The series of peaks at high f2 values all correspond to various 
combinations of bond angles in the (disordered) triangular phase. The interfacial region has corresponding to the first dip 
in the curve, (b) Illustration of the spatial distribution of f2 values as stated in (a) f2 < 0.2 (square) are black, those for 
1.3 > n > 0.2 are red (interface) and those for Q. > 1.3 are blue (triangle), (c) The fluctuation of the local velocity about 
the centre of mass value < (?; — Vcm)^ > plotted against Q. It is clear that the fluctuations are largest for the particles in the 
interfacial region, (d) The fluctuation of the local velocity as a function of the local volume per particle. 



direction are always long-ranged while those along the 
transverse direction are generically short-ranged. Fifth, 
the intuitively unexpected expansion of the plastic flow 
regime at large and its subsequent collapse. Finally, we 
have suggested a possible connection to the classic prob- 
lem of the origin of peak effect anomalies, conjecturing 
that the explanation might be found in the very nature 
of the unusual dynamical states obtained when a system 
close to a first-order (structural or melting) phase transi- 
tion is driven across a quenched-disordered background. 

There are several realizations of non-equilbrium states 
of complex fluids, such as sheared lamellar phases , worm- 
like micelles, and driven colloidal suspensions, which ex- 
hibit remarkably non-trivial behaviour as a consequence 
of dynamical phase transitions [8l|, [H, [sl]. Possible rela- 
tions to equilibrium descriptions of systems at phase co- 
existence have also been outlined[84]. Many features, in- 
cluding profoundly non-linear response, complex spatio- 
temporal behaviour and large noise signals appear to 



be common to these systems when driven away from 
equilibrium[81] . What singles out the category of systems 
we study here is the additional complication of quenched 
disorder. Clearly these are suggestive links at the in- 
tersection of these various fields. Further exploration of 
these connections would appear to be fruitful. 
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FIG. 19: [Color Online:] Configuration snapshots arising from a quench to zero force of a randomly chosen configuration within 
the coexistence regime, at time steps (a) 0, (b) 4x10^, (c) lO'', (d) 5 xlO'', (e) 10^ (f) 2 x 10^, (g) 3 x 10^, (h) 5 x 10^ and 
(i) 9 X 10^, with vs = 6.0. The system size is N = 10000. 
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